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We present a new tool for relating theory and experiment suited for non-Gaussian 

,.^J ' theories: non-Gaussian spectra. It does for non-Gaussian theories what the angular 

QJ ' power spectrum Ci does for Gaussian theories. We then show how previous studies 

^^ ' of cosmic strings have over rated their non-Gaussian signature. More realistic 

maps are not visually stringy. However non-Gaussian spectra will accuse their 
stringiness. We finally summarise the steps of an undergoing experimental project 
aiming at searching for cosmic strings by means of this technique. 

1 Theory and experiment in CMB physics 

In this review we present a new tool for making contact between theory and 
CNJ ' experiment in cosmic microwave back;ground (CMB) physics Id. The central 

^^ , tool which has played this role in the past is the angular power spectrum 

(known as Ci). If the underlying theory is Gaussian it is known that the 

power spectrum C^ encodes the totality of the predictions made by the theory. 

i-Q ' It then makes sense to direct all experimental effort towards its measurement. 

Also the power spectrum represents the data-reduction end-point should we 

O ' firmly believe in the Gaussianity of the underlying Universe. If we start with 

fi , an all sky map with 10^ pixels one could reduce them to roughly lO'^ values 

"^ ' for Ci. For Gaussian theories all the science is in these 10"^ numbers. The map 

itself is redundant. If we further believe not only that the Universe is Gaussian, 

but also that inflation is the truth, then these Cp are only dependent on about 

10 parameters. Therefore we may in fact carry this reduction further to about 

;^ ' 10 quantities. 

If the underlying theory is non-Gaussian, however, the power spectrum is 
not the end of the story. In this paper we define a new set of spectra, which we 
will label non-Gaussian spectra, and which carry predictive power in generic 
non-Gaussian theories. We will argue that non-Gaussian spectra could in fact 
be the best arena for connecting theory and experiment in motivated non- 
Gaussian theories such as cosmic string scenarios. Although the estimation of 
these spectra complicates data-reduction and data-analysis considerably, non- 
Gaussian spectra are essential for bringing out the full predictive power of non- 
Gaussian theories. They may also help to quantify in which sense topological 
defect theories are non-Gaussian. This becomes particularly relevant if one 
allows for the full complication of theories like cosmic strings to be included in 



the consideration of their non-Gaussian signals. 

It should be stated from the start that if the underlying theory is indeed 
Gaussian, non-Gaussian spectra are perfectly useless. One may not, for in- 
stance, improve estimates of cosmological parameters in inflationary scenarios 
through their measurement. Prejudice is however a nasty element in science, 
and one should think of non-Gaussianity as an issue in its own right, regard- 
less of any theoretical dogma. Most data-analysis methods, in particular, rely 
blindly on Gaussianity. They could well break down miserably should the data 
prove to be non-Gaussian in the first place. 

2 A taste of non-Gaussian spectra 

The idea of non-Gaussian spectra is to supplement the angular power spectrum 
with a set of quantities whose number may never exceed the number of infor- 
mation degrees of freedom one starts from. If N independent pixels have been 
measured then transforming them into a larger number of quantities obviously 
introduces redundancy. On the other hand computing a small number of quan- 
tities, such as the skewness and kurtosis, or in general submitting the data to 
Gaussianity tests, may leave room for all sorts of perverse and less perverse 
non-Gaussianity to pass unnoticed. In general our philosophy is to factor out 
from the data only rotational degrees of freedom (how we oriented the axes) . 
There arc 3 degrees of freedom in rotations. There should then be A^ — 3 
quantities with which to compute the power spectrum, plus the non-Gaussian 
spectra. 

Let us consider the spherical harmonic expansion of an all-sky map: 

1=0 m=-£ 

Consider also the situation where the field is very small so that instead one 
may expand in Fourier modes: 



AT(x) 
T 



^J^aiky^- (2) 



There is some correspondence between af^ and a(k). To some extent the 
moduli |k| correspond to £, and act as the scale of the modes. On the other 
hand the to is a bit like the direction of the vector k and labels something 
like the direction of the mode. One must remember that these modes are 
complex, and that therefore have phases as well as moduli (which are the only 
components entering measures of power) . 




Figure 1: A schematic picture of power spectra and non-Gaussian spectra. Dividing Fourier 
space into rings with AA; = 1 one obtains the power spectrum Ci by averaging the power 
on each ring. Non-Gaussian spectra Bm are a measure of how the power is distributed in 
angle, in each ring. Phase spectra and inter-ring spectra complement this, so-called, shape 

spectra. 



If for simplicity one concentrates on the Fourier space picture, then com- 
puting the power spectrum may be seen as the result of the following algorithm. 
One divides Fourier space into rings with Afc = 1, and then averages the square 
of the moduli of modes lying in each of these rings. The average Ce (or C{k)) 
represents a measure of the amplitude (or the power) of the fluctuations on 
the scale i (or |k|): 

C(fc) = (|a(k)|2) (3) 

Within this picture one can think of a natural way of extending this construc- 
tion as a set of transverse or ring spectra. These represent a measure of how 
the power is distributed across each ring in angle. Such a measure of angular 
distribution of power would naturally represent the shape of the fluctuations. 
However one may now talk about shape on a given scale, something which is 
far from visually intuitive. 

On top of this one can consider phase spectra. Phases transform under 
translations and therefore represent the localisation of the fluctuations. Again 
we may define localisation as a function of scale, something rather abstract. 
Finally one can consider a radial spectrum of correlations between adjacent 



rings. Intcr-ring correlators are the crucial aspect of non-Gaussianity which 
allows for shape and localisation on the various scales to be transmuted into 
structures which we can recognise visually. Indeed one needs the constructive 
interference between all the scales so that something as abstract as shape and 
place on a given scale becomes, say, the picture of an elephant. 

The advantage of such abstract definitions of shape and place is precisely 
in that they do not rely on visual recognition. As we will see soon very rarely 
life is so kind to us so as to provide us with evident visual non-Gaussianity. A 
more abstract, but still comprehensive, method for describing shape and place 
is therefore necessary. 

3 Realistic stringy skies 

The formalism we have just outlined acts as comprehensive formalism for en- 
coding non-Gaussianity, which is an advantage over Gaussianity tests. Also, 
unlike the n-point correlation function, it does not contain redundant infor- 
mation. However it lives naturally in Fourier space. This contradicts a well 
established dogma: non-Gaussianity is obvious in real space, but gets diluted 
in Fourier space. According to this argument a Fourier mode is the result of 
the addition of many possibly non-Gaussian temperatures: 

a(k)=/|^e-- (4) 

The central limit theorem then tells us that even if the — J^ are very non- 
Gaussian the a(k) will be approximately Gaussian. As an example we may 
consider the Kaiser-Stebbins effect from cosmic strings (see Fig. 0, top figure). 
Distinctive stringy discontinuities may be recognised in the top map. The 
Fourier transform of such a map would be a bit of a mess, a mess that is not 
particularly non- Gaussian. 

We will counter this argument using precisely the example of cosmic string 
maps. These have been oversimplified in the past, in a way that over rates their 
non-Gaussinity. A more realistic examination of stringy skies does not comply 
with the cubist microwave sky often attributed to cosmic strings. To begin 
with realistic cosmic strings are rather contorted objects. As the top figure 
in Fig. @ shows one does see jumps on stringy skies but these have a rather 
limited coherence length. 

More important than this is the recognition that even the top picture in 
Fig. |2| is a simplification. It assumes that the glow of photons coming out 
of the last scattering surface is perfectly homogeneous. However one must 
remember that there were strings before last scattering. These also caused 




Figure 2: The Kaiser Stebbins effect from cosmic strings (top) if the glow out of the last 

scattering surface were homogeneous. Adding on the perturbations induced by strings before 

last scattering leaves what looks like a Gaussian map (bottom). 
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Figure 3: Schematic representation of the power spectrum of the Gaussian fluctuations 
induced by cosmic strings before last scattering and the non-Gaussian fluctuations induced 
after last scattering by the Kaiser-Stebbins effect. We also indicate the scale of the inter- 
string separation at last scattering, which is 9 ^ 20' corresponding to ^ R:i 60° / 6 ^ 200. 



fluctuations in the cosmic radiation. After last scattering radiation is made 
up of free photons that fly past the moving strings and "take a picture" of the 
string network. Before last scattering the photons behaved more like a fluid, 
which was slashed by the strings, but was also subject to its own pressure. A 
complete mess of waves may be expected, rather than a neat picture of the 
string network. Indeed the only available calculation of string pierturbations 
before last scattering shows rather Gaussian looking fluctuations □. 

One may plot the power spectrum of the fluctuations induced by cosmic 
strings schematically as in Fig. |[ The Ci may be seen as the result of two 
contributing components. One corresponds to fluctuations induced by strings 
after last scattering, which gives the beautiful Kaiser-Stcbbins effect one usu- 
ally pictures. The other corresponds to the fluctuations induced before last 
scattering, which should be roughly Gaussian. The non-Gaussian component 
dominates on large scales, presenting a slightly tilted spectrumQ, then falls off 
as a power lawQ. The Gaussian component is negligible at large scales (white 
noise, rather than quasi-scale invariant), then rises into a singla-Doppler peak 
of height as yet unknown, but roughly placed at £ ~ 400 — 600 Q. It then falls 
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Figure 4: The Fourier space structure of realistic string skies. 



off exponentially, due to Silk damping ^. 

If the Gaussian component were not there all one would need would be 
an experiment with a resolution better than the inter-string separation at last 
scattering, and the Kaiser-Stebbins effect would emerge in all its glory. How- 
ever at the inter-string separation scale the power of the fluctuations is domi- 
nated by the Gaussian component. Still not all is lost. Because the Gaussian 
component falls off as an exponential, whereas the non-Gaussian component 
falls off like a power law, at very high i the fluctuations should be dominated 
by the non-Gaussian component again. All in all the Fourier space image of 
a stringy sky may be divided into mainly Gaussian and mainly non-Gaussian 
bands, as depicted in Fig. ^. It looks as if realistic string maps lead to the 
conclusion that one needs much higher resolution to detect strings than previ- 
ously thought, and that their detection will probably be clearer in the Fourier 
domain. 

We are now in position to reverse the central limit theorem argument so 

"We believe that the point made inEl in purely semantic, concerning what to call after and 
before last scattering. Whatever the before/after definition is, the power in the Gaussian 
component should always fall off exponentially 



often used against non- Gaussian! ty analysis in the Fourier domain. The real 
space temperature is now seen as the sum of many Fourier space components, 
some of which are strongly non-Gaussian, and some others nearly Gaussian: 

AT(x) _ f rfk^^j^^^,,., ^^^ 



T J 2tt 

The central limit theorem can now, even more than before, predict that the 
real space imagine will be nearly Gaussian. The bottom figure in Fig. shows 
that indeed the realistic Kaiser-Stebbins effect does not lead to any structures 
which any conceivable Human could visually recognise. 

Clearly we have a situation where the Gaussian and non-Gaussian com- 
ponents are naturally separated in Fourier space. The non-Gaussian signal in 
Fourier space will never be a visual signal. However it should be now clear that 
if I find an algebraic way to detect stringiness in Fourier space, this method 
will be robust against the addition of Gaussian signal or even Gaussian noise, 
as long as a non-Gaussian band in Fourier space survives these additions. 

Non-Gaussian spectra provide such an algebraic tool. They define shape on 
a given scale. They should therefore allow the detection of stringy shapes on the 
outer non-Gaussian band where the stringy Kaiser-Stebbins effect dominates 
the signal. 

4 Non-Gaussian spectra 

A more mathematical definition of non-Gaussian spectra will now be given. 
Consider a ring of the Fourier space where Nk independent complex modes 
a(ki) = 3?[a(ki)] -|-i3[a(ki)] live. In Gaussian theories these are distributed as 

F(5R[a(kO],3Kk.)]) = ^^^^^i,^x 

^^P- [^E(K'[«(k.)] + 3'[a(k.)])j (6) 

where mk — Nk/2. First separate the Nk complex modes into ruk moduli pi 
and TTifc phases <pi 

K[a(ki)] = Pi cos (j>i 

3[a(ki)] = Pisin^i (7) 

The Jacobian of this transformation is 



a(K[a(k,)],3[a(k,)]) 



Cl{piAi) 



1=1 



The {pi} may be seen as Cartesian coordinates which we transform into polar 
coordinates. These consist of a radius r plus rrik — 1 angles 9i given by 

Pi = r cos 9i TT sin 9j (9) 

3=0 

with sin 6*0 = cos^m^ = 1. In terms of these variables the radius is related to 
the angular power spectrum by C{k) — r^ /{2mk)- In general the first mk — 2 
angles 9i vary between and tt and the last angle varies between and 2tt. 
However because all pi are positive all angles are in (0, 7r/2). The Jacobian of 
this transformation is 

Q/ N nife-l 

0[Pl,---,PmJ 



d{r, 



1 H sin™''-*0Vi (10) 



4 = 2 



Polar coordinates in ruk dimensions may be understood as the iteration of the 
following rule: 

Pi = n cos 9i 
r,_i = riSin6'i (11) 

in which r^ is the radius of the shade mk — i + I dimensional sphere obtained 
by keeping fixed all pj for j = 1, • • • , i — 1: 



p1 + p1+i + --- + pI, (12) 

One may easily see that this is how 3D polars work, and also that the transform 
(^ follows this rule. Hence one may invert the transform (g) with 

9i — arccos — (13) 

V pf + p?+i + ■ ■ • + pIi, 

for i = 1 , • • • , TOfc — 1 . 

The Jacobian of the transformation from (5R[a(ki)], S[a(ki)]) to {r, 9i, (j)i} 
is just the product of (0) and iRw. Hence for a Gaussian theory one has the 
distribution 

r^'=-iexp-ff^) "'-c-i 
n^'^^'^^) - (2^,2)^. V n cos0~.(sin0^)^---i (14) 
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In order to define 6i variables which are uniformly distributed in Gaussian 
theories one may finally perform the transformation on each 6i: 

9, ^ sm^'-^\e,) (15) 

so that for Gaussian theories one has: 

^ ' i—l 

The factorization chosen shows that all new variables are independent random 
variables for Gaussian theories, r^ has a Xn distribution, the "shape" vari- 
ables 0i are uniformly distributed in (0,1), and the phases 4>i are uniformly 
distributed in (0, 27r). 

The variables 9i define a non-Gaussian shape spectrum, the ring spectrum. 
They may be computed from ring moduli pi simply by 

They describe how shapeful the perturbations are. If the perturbations are 
stringy then the maximal moduli will be much larger than the minimal mod- 
uli. If the perturbations are circular, then all moduli will be roughly the same. 
This favours some combinations of angles, which are otherwise uniformly dis- 
tributed. In general any shapeful picture defines a line on the ring spectrum 9i . 
A non-Gaussian theory ought to define a set of probable smooth ring spectra 
peaking along a ridge of typical shapes. 

We can now construct an invariant for each adjacent pair of rings, solely 
out of the moduli. If we order the pi for each ring, we can identify the maximum 
moduli. Each of these moduli will have a specific direction in Fourier space; 
let hmax and k^^^^ be the directions where the maximal moduli are achieved. 
The angle 

■0(fc, k') = -a,ng{k^ax,'k'^ax) (18) 

TT 

will then produce an inter-ring correlator for the moduli, the inter-ring spectra. 
This is uniformly distributed in Gaussian theories in (—1,1). It gives us in- 
formation on how connected the distribution of power is between the different 
scales. 

We have therefore defined a transformation from the original modes into 
a set of variables {r, 0, (/>, i/'}- The non-Gaussian spectra thus defined have a 
particularly simple distribution for Gaussian theories. The fact that this dis- 
tribution does not have a peak shows clearly that we cannot use non-Gaussian 
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spectra to find out anything about tlic parameters of a Gaussian theory (eg. 
cosmological parameters in inflation). We shall call perturbations for which the 
phases are not uniformly distributed localized perturbations. This is because 
if perturbations are made up of lumps statistically distributed but with well 
defined positions then the phases will appear highly correlated. We shall call 
perturbations for which the ring spectra are not uniformly distributed shapeful 
perturbations. We will identify later the combinations of angles which measure 
stringy or spherical shape of the perturbations. This distinction is interesting 
as it is in principle possible for fluctuations to be localized but shapeless, 
or more surprisingly, to be shapeful but not localized. Finally we shall call 
perturbations for which the inter-ring spectra are not uniformly distributed, 
connected perturbations. This turns out to be one of the key features of stringy 
perturbations. These three definitions allow us to consider structure in various 
layers. White noise is the most structureless type of perturbation. Gaussian 
fluctuations allow for modulation, that is a non trivial power spectrum C{k), 
but their structure stops there. Shape, localization, and connectedness con- 
stitute the three next levels of structure one might add on. Standard visual 
structure is contained within these definitions, but they allow for more abstract 
levels of structure. 



5 The meaning of non- Gaussian spectra 

The decomposition {C(fc),0, 0, ^} has an immediate physical interpretation. 
The angles 9 reflect the angular distribution of power, and therefore reflect 
shape. The phases (j) transform under translations and so contain the infor- 
mation on position and localization of the structures in the field. The angles 
■0 correlate different scales, and therefore tell us how connected the structures 
are. For a Gaussian random field the variables {9, (p, V'} are all uniformly dis- 
tributed reflecting complete lack of structure besides the power spectrum. In 
terms of the various levels of structure considered we can then characterize 
Gaussian fluctuations as shapeless, delocalized and disconnected. By compar- 
ison with a Gaussian we may then define structure at different levels. We will 
say that fluctuations for which 9 are not uniformly distributed are shapeful. 
If the (j) are not uniformly distributed we shall say the fluctuations are local- 
ized. If the V' are not uniformly distributed the fluctuations are connected. 
Although visual structure has room within these definitions, they are consid- 
erably more abstract and general. We may consider highly non visual types 
of structure such as shapeful but delocalized fluctuations or disconnected lo- 
calized stringy fluctuations. In this sense we regard our formalism as a robust 
definition of structure, which goes beyond what is visually recognazible and so 
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Figure 5: A spherical hot spot which has been deconstructed at different levels. On the 
top left hand panel we have the pure non-Gaussian signal. The angles 8i have been redrawn 
uniformly on the top right picture. On the bottom left the phases (pi were redrawn unformly. 
On the bottom right we applied an independent unformly distributed rotation on all rings 
in Fourier space. From top to bottom and left to right, a plain regular sphere, a shapeless 
sphere, a delocalized sphere, and a disconnected sphere. 
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Figure 6: The Kaiser-Stebbins effect (top left) and its various stages of deconstruction. The 
angles 6i have been redrawn uniformly on the top right picture. On the bottom left the 
phases </<; were redrawn unformly. On the bottom right we applied an independent unformly 
distributed rotation on all rings in Fourier space. Respectively we have strings, shapeless 
strings, unlocalized strings and disconnected strings. 



13 



is tied down to our particular and narrow path of natural selection. We may 
imagine an alien civilization with Fourier space eyes (say interferometric eyes 
Q), and a brain trained to recognize Fourier space structure at many different 
levels, structure that would seem totally non obvious to our human eyes. 

To illustrate the limitations of human vision we shall now destroy highly 
structured maps level by level, that is Gaussianize only one of the variable types 
{9,(f>,ip}. Initially there will be structure at every level, shape, position, and 
connectedness. We will remove structure gradually, a fact not disasterous for 
the alien civilization referred above, but which will illustrate the limitations 
of the human visual method for recognizing non-Gaussianity. In Figure H 
we play this game with a sphere. We depict a spherical hot spot in real 
space, then a shapeless sphere, a delocalized sphere, and a disconnected sphere. 
For the case of a sphere we find that what we recognize as shape is mostly 
localization. A shapeless sphere keeps its recognizable features. On the other 
hand a delocalized sphere loses it characteristic features. Indeed the idea of 
a shapeful but non-localized object sounds somewhat surreal for all we can 
visually conceptualize. Nevertheless our formalism will acuse the strong but 
not obvious non-Gaussianity exhibited by a delocalized sphere. 

In Figure we repeat the same exercise for a map displaying the Kaiser- 
Stebbins effect from cosmic strings. Shapeless strings, delocalized strings, and 
disconnected strings are shown. Considerable disarray is introduced in every 
case, but one may say that disconnected strings as well as delocalized strings 
are perhaps the most messy of them. This is consistent with the strong signal in 
-0 we have found for the case of the realistic Kaiser Stebbins effect. On the other 
hand the fact that line- like discontinuities are present even for shapeless strings 
shows how much more structure there is in the map on top of the structure 
which we can recognize. This is important since the beautiful patchwork is 
very fragile to the hard realities of noise and supperposed Gaussian signal. In 
the real world, it turns out, the non-visual feature which is the connectedness 
of strings happens to survive much better than the patchwork (which reflects 
mostly localization). 

6 Non-Gaussian spectra and cosmic string detection 

As shown in Sec. 3 the existence of a Gaussian background in string scenarios 
implies the need for higher resolution than thought before in order to de- 
tect cosmic string non-Gaussianity. Interferometer experiments are therefore 
favoured. For such high resolution interferometers can only look into small 
fields. We therefore concentrate on a field of some 20' across. The Kaiser- 
Stebbins effect before and after the addition of the Gaussian signal would then 
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Figure 7: A field of the order of 20' across sliowing tfie Kaiser-Stebbins effect from a cosmic 
string before and after the Gaussian signal is added. 

look like Fig. |. 

Clearly the string is no longer visible. What is more, we have checked that 
the traditional non-Gaussianity tests fail to detect the hidden small scale non- 
Gaussianity in the right picture in Fig. 0. We have checked this statement with: 
pixel histograms, skewness and kurtosis, density of peaks, topological tests, 
the 3-point correlation function, and other traditional tests. As an example we 
show in Fig. g the failure of the Euler characteristic test to detect the realistic 
Kaiser-Stebbins effect. 

However we found that non-Gaussian spectra may still detect the hidden 
string after Gaussian signal has been added. Let us first look at shape spectra 
(Fig. ^. For rings at low k these are very Gaussian. However, as we go out 
into high k rings, something peculiar starts to happen. The straight string 
spectrum is of course never visible. However a clear ridge in the distribution 
of shapes emerges. This corresponds to the wiggly string shape. There is 
naturally cosmic/sample variance in the shape spectrum. This is because not 
all wiggly strings have the same wiggles. Therefore, as in the case of Ci spectra, 
one must complement the average spectrum with error bars. Still the concept 
of shape spectrum makes sense, because shapes are not uniformly distributed, 
but rather have a peaked distribution. 

More impressive still is the inter-ring spectrum, shown in Fig. nO. If an 
object is exactly stringy then it concentrates the power in Fourier space along 
one direction. More importantly, this direction is the same for adjacent rings. 
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Figure 8: The mean Euler characteristic, F as a function of threshold for a string map with 
Gaussian noise (sohd line) and for a pure Gaussian map with the same power spectrum (the 
shaded region is the 1 a region around the Gaussian mean, estimated from a 100 realizations). 
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Figure 9: The ring spectrum at high k inferred from 100 realizations. The shaded region 

represents a probability larger than - for the the values of 9i to occur (this generalises the 

concept of 1 sigma error bar for non-Gaussian distributions. 
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Figure 10: The inter-ring spectrum ^. The shaded region represents a probabiUty larger 
than - for the the values of ipi to occur. 



Then, for wiggly strings, one may expect that the direction where the maxi- 
mum modulus is achieved is highly correlated from ring to ring. As we see in 
Fig. |l^, as we go up in k the angle ijj between maximal directions is at first 
uniformly distributed. Gradually, as the string signal starts to dominate, the 
-0 distribution becomes highly peaked around ^ = 0. This is a very impressive 
signal, with a very small error bar. We put most of our hopes of detection on 
this signal. 



7 Prospects for the future 

The statistics presented here should be comprehensive in their detection of 
non-Gaussianity, but naturally we are concentrating on finding characteristic 
signatures for string networks in this framework. 



7.1 Simulating the sky 

We plan to test our statistics on actual observations. But before using real 
data, we need to test the method on simulations. We are currently working 
on simulated brownian strings. These are useful in testing our intuition about 
the statistics but of course are not the most realistic examples of strings. The 
most accurate string maps are generated by evolving a string network from its 
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initial conditions up to the present. We will do this by means of the integer 
string code described inli3. We hope to generate 10000 maps in this way. 

Once we have a stri^jg map, the photon temperature map is calculated 
using Laplace's equations 

AT ^ ^ 

V^— = -87rGV-C/ (19) 

where all operations here arc on the plane orthogonal to some direction vector 
n. If n is chosen to be in the z direction then we have 

/ + 00 
[Qo,{t,X^{t)) - e,3{t,X^{t))]dt (20) 

-oo 

where O^i^ is the stress-energy tensor for a string and i — 1,2. Laplace's 
equation is most easily solved in Fourier space so we are like an interferometer 
experiment, which observes directly in Fourier space. 

We will assume the fluctuations on the last scattering surface to be Gaus- 
sian, so all we need to specify is the power spectrum. As the height of the peak 
in the spectrum is unknown for defect models and its position is uncertain we 
need to fit parameters to the C^. 



7.2 Simulating the experiment 

First we have to take into account effects of limited sky coverage and reso- 
lution. This is done for an interforemeter by modelling the primary beam of 
the antennae with a window function W(x.) in real space and the synthesized 
beam by a function B{x). The Fourier transform of B{x), i?(k), is 1 where 
the Fourier domain is sampled and otherwise so it is effectively a window in 
Fourier space. Therefore the limit of resolution in real space places an upper 
limit on the values of k we can look at. Large values of k are precisely what 
we want to look at, so we hope for a high resolution experiment. The Fourier 
transform of W(x), WCk), corresponds to a limit in the resolution in Fourier 
space since Vt^(x) cuts off large scales in real space. This introduces correla- 
tions in the a(k), so we must be careful in deciding how to sample Fourier 
space. 

Foreground emission can be split into two main components: smooth emis- 
sion from dust, synchrotron and free-free emission and point sources. Smooth 
emission gives approximately constant errors for all k so it can be modelled as 
Gaussian noise (see below) . Spherical point sources have pi = constant in each 
ring, so they should not interfere with the ring spectra of strings. 
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Finally we need to add noise to our sky map. This is assumed to be 
Gaussian. It adds an extra term to the covariance matrix of the Fourier modes 
a(k) which may be written asQ 

(a^(k.),a^(k,))=5.,^^^^^ (21) 

with uj~^ = j^-, where A is the area of Fourier space covered, s^ is the 
sensitivity of the detector, itot is the total observing time, Uf is the number of 
fields observed and fl is the area of each field. We will experiment with the 
various parameters to find the best observing strategy. 

7.3 The statistics 

Once we have our maps of stringy skies we can test the statistics on them so 
we will be in a position to decide what the best observing strategies are {ie 
how long should be spent in each field, how many fields should be observed 
etc). 

Other possibilities for statistics to investigate include: 

• Quantities made up of the phases of the a(k) (in the same way that the 
shape spectra are made up of the moduli). The phases contain all the 
information about position in real space and so are complementary to 
the moduli. 

• Looking at the two point function around a ring in Fourier space. 

• Taking another Fourier transform in angle around a ring in Fourier space. 
This and the two point function are alternative ways of looking at the 
distribution of power in rings. 

Eventually we hope to test our statistics on data from the VLA and BIMA 
experiments. These are interferometer experiment with a field of a few arcmin- 
utes but very high resolution and sensitivity. 
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